Coexistence of distinct charge fluctuations in 6*-(BEDT-TTF) 2 X 
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Using the Lanczos exact-diagonalization and density-matrix renormalization group methods, we 
study the extended Hubbard model at quarter filling defined on the anisotropic triangular lattice. 
We focus on charge ordering (CO) phenomena induced by onsite and intersite Coulomb interactions. 
We determine the ground-state phase diagram including three CO phases, i.e., diagonal, vertical, 
and three-fold CO phases, based on the calculated results of the hole density and double occupancy. 
We also calculate the dynamical density-density correlation functions and find possible coexistence 
of the diagonal and three- fold charge fluctuations in a certain parameter region where the onsite and 
intersite interactions compete. Furthermore, the characteristic features of the optical conductivity 
for each CO phase are discussed. 

PACS numbers: 71.10.Fd, 71.20.Rv, 71.30.+h, 71.45.Lr 
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I. INTRODUCTION 

In a number of correlated electron systems, ordering of 
charge degrees of freedom of electrons, or charge order- 
ing (CO), is a key feature to understanding their physi- 
cal properties because it is intimately related to the elec- 
tronic transport and magnetic phenomenal The concept 
of CO was originally introduced to interpret an array or- 
der of Fe 2+ and Fe 3+ ions in magnetite Fe304, which 
is known as the Verwey transition^ For several years, 
the importance of CO has been increasingly recognized 
in conjunction with spatial localization of charge carriers 
observed in manganese oxide Lai_ K Ca x Mn03 (x > 0.5) 
copper oxide Lai.6-xNdo.4Sr a: Cu04 (x = 0.12)^ etc. Es- 
pecially, in the field of low-dimensional organic materi- 
als, the CO itself has been one of the main issues for 
discussion^ Quite recently, peculiar charge fluctuations 
have been found in the CO state of quasi-two-dimensional 
(2D) organic conductors 6>-(BEDT-TTF) 2 X (hereafter, 
BEDT-TTF is abbreviated as ET) and they have at- 
tracted much attention due to their intriguing relation 
to the superconducting (SC) state. 

The crystal structure of #-(ET) 2 X are built with 
an alternating stack of ET conducting and anion (X) 
insulating layers. In the conducting layer, the ET 
molecules form an anisotropic triangular-lattice structure 
(see Fig. [TJ and the band is |-filled in terms of electrons 
or ^--filled in terms of holes. In experiments, the transfer 
integrals (t p ,t c ) are controllable by altering a dihedral 
angle (6) with the substitution of anion and/or adjust- 
ment of pressure.^ A decrease in the dihedral angle cor- 
responds to an increase in t p and therefore a decrease in 
t c /t p . At low temperatures (T < 250 K), the compounds 
exhibit a variety of phases, such as CO insulator, para- 
magnetic metal, and SC, function of the dihedral 
angle. Overall, the metal-insulator transition tempera- 
ture (Too) decreases with decreasing the dihedral an- 
gle. Only 6*-(ET) 2 l3 that has the smallest dihedral angle 
9 = 100° in the 0-ET family shows SC in place of CO 
below T w 3.6 K4 



Thus far, the CO pattern in each compound has 
been investigated by X-ray diffraction and 13 C-NMR 
measurements. In the series of X= MM' (SCN)4 with 
M' =Co and Zn, M =T1 and Rb, a stripe-type CO 
which has a twofold periodicity along the c axis is ob- 
served below Tco^^^- The compounds with M =T1 
and Rb is known to have an average dihedral angle of 
6 = 120° to 110°. Using the estimated ratio of the 
transfer integrals t c /t p ~ 0.5 — 0.4, the stripe-type CO 
was successfully reproduced by the mean-field theory^ 
But then, the X-ray diffuse scattering experiment a 14 : 15 
demonstrated the presence of peculiar charge fluctua- 
tions in 6>-(ET)CsM'(SCN) 4 with M =Co and Zn; two 
diffuse peaks with the wave numbers qi = (q a ,qi,,q c ) = 
(0, k, 1/2) and q 2 = (2/3, k, 1/3) were found. The wave 
number q! is the same as that of the stripe-type CO ob- 
served in the compounds with M =Tl and Rb, but the 
other wave number q 2 corresponds to a three-fold period- 
icity. The compounds with M =Cs have relatively small 
dihedral angle 8 = 104° and are located in the vicinity 
of the quantum critical point Tco ~ K. Accordingly, 
the ratio of the transfer integrals is expected to be rather 
small, i.e., t c /t p <C 1, which is different from the case 
of M =T1 and Rb. Thus, a very absorbing subject, i.e., 
the coexistence of charge fluctuations with different wave 
numbers, has been provided. 

The appearance of the spatially inhomogencous CO 
has also been suggested by the dielectric response, ac re- 
sistivity^ 1 H-NMR, EPR, and static magnetic suscepti- 
bilit y 17 ' 18 measurements. Furthermore, a power-law be- 
havior of the current-voltage characteristics over a wide 
range of currents^ and a current-induced melting of in- 
sulating CO domains 1 ^ have been observed. In addition, 
the optical conductivity spectra show a transition from 
short-ranged to long-range CO around T = Tco, an d 
stronger insulating features appear in MM' =CsZn than 
in MM' =RbZni 21 i 22 A number of theoretical studies 
have then been carried out in order to clarify the nature 
of the CO and charge fluctuations^ 3 .^^^^^ 

Motivated by such developments in the field, we con- 
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nearest-neighbor pairs, tij and Vij are hopping integral 
and intersite Coulomb repulsion between sites i and j; 
here we retain two kinds of nearest-neighbor hopping in- 
tegrals (repulsions) t p (V p ) and t c (V c ) as shown in Fig.Q] 
U is the onsite Coulomb repulsion. We restrict ourselves 
to the case at quarter filling; i.e., (rii) — 0.5 where (• • •) 
denotes the ground-state expectation value. In this pa- 
per, we consider the case of U/t p = 10 and t c /t p = as a 
typical set of parameter for X=CsCo(SCN)4^ Hereafter, 
we take t p = 1 as the unit of energy. 



FIG. 1: (Color online) Schematic representation of the 
anisotropic triangular-lattice Hubbard model. Each circle de- 
notes an ET molecule. 

sider in this paper the extended Hubbard model at quar- 
ter filling defined on the anisotropic triangular lattices. 
We employ the density-matrix renormalization group 
(DMRG) and Lanczos exact-diagonalization methods to 
investigate the electronic states of the model. First, we 
calculate the hole density and double occupancy to deter- 
mine the ground-state phase diagram. Next, the dynam- 
ical density-density correlation functions are calculated 
to study the low-energy charge excitations. We also cal- 
culate the single-particle excitation spectrum to elabo- 
rate the anomalous metallic states seen in the low-energy 
charge excitations of the three-fold CO phase. We finally 
obtain the optical conductivity for each phase and clar- 
ify the behavior of the charge degrees of freedom of the 
model by focussing on the three-fold charge fluctuations. 

This paper is organized as follows. In Sec. II, the 
2D extended Hubbard model on the anisotropic triangu- 
lar lattices is introduced. We also define some physical 
quantities and explain the applied methods for the cal- 
culations. In Sec. Ill, we present the ground-state phase 
diagram of the extended Hubbard model and discuss the 
physical properties related to the CO in each phase. Sec- 
tion IV contains conclusions and discussions including 
comparison with the experimental results. 



II. MODEL AND METHOD 
A. Hamiltonian 



B. Physical quantities 

We employ two kinds of numerical methods, i.e., the 
DMRG and Lanczos methods. Either of these methods is 
chosen for the calculation of each physical quantity. The 
details are explained in the following. 

1. Static quantities 

We are interested in the CO phenomena in the subject 
materials, so that it is useful, first of all, to make the 
ground-state phase diagram associated with the charge 
degrees of freedom. To this end, we calculate the hole 
density (rii) and double occupancy (di) (di = n^n^) for 
all sites i to investigate the charge distribution. A large 
size system is required for reproducing all possible CO 
patterns and for an accurate estimation of these quanti- 
ties. We thus study a finite-size cluster of L a = 8 and 
L c = 6 using the DMRG method. The periodic boundary 
conditions (PBC) are applied for the c-direction, whereas 
the open-end boundary conditions (OBC) are applied for 
the a-direction. This choice of the boundary conditions 
enables us to detect the CO state clearly (see Sec. Ill A). 
We keep up to m ~ 1400 density-matrix eigenstates in 
the DMRG procedure and thus the maximum truncation 
error, i.e., the discarded weight, is less than 1 x 10~ 4 . In 
this way, the maximum error in the ground-state energy 
is estimated to be AEo/N ~ 10~ 2 for the most inaccu- 
rate case V p = V c = 0. Based on the results of the hole 
density and double occupancy, we determine the ground- 
state phase diagram on the parameter space (V p , V c ). 



The Hamiltonian of the extended Hubbard model de- 
fined on the anisotropic triangular lattice is given by 

H= - 5^ty(4^+H.c.) 

i (v) 

where cj CT {ci a ) is the creation (annihilation) operator of 
a hole with spin a at site i, ni a — cl a Ci a is the number 
operator, and n; = + n^. The sum (ij) runs over 



2. Dynamical density- density correlation 

Next, we consider the dynamical density-density cor- 
relation function defined as 

AT(q, W ) = Y, I (€ h \n^ h )\ 2 S(u; - E v + E ) (2) 

in order to study the low-energy charge excitations. Here, 
\tp^ h ) and E^ h are the f-th eingenstate and eigenen- 
ergy of the system with Nh holes iv = corresponds 
to the ground state). Since the exact definition of the 
momentum-dependent operators with the OBC is quite 
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difficult, we choose the PBC for a quantitative estimation 
of the spectrum; we therefore use the PBC for both the x- 
and y-directions (or the a- and c-directions). With this 
PBC, the density operators n q can be precisely defined 
by 

v i 

where N is the number of lattice sites L x x L y (or L a x 
L c ) and is the position of site i. However, in this 
case, it is quite difficult to carry out sufficiently accurate 
calculations with the DMRG method, so that we here use 
the Lanczos method on small clusters with 16 (L x = 4, 
L y = 4) and 18 (L a — 6, L c = 3) sites. Consequently, 
the system size is restricted but the result is numerically 
exact. 



3. Single-particle excitation spectrum 

Of interest are also the evolution of the fundamental 
excitations in the systems with strongly frustrated cor- 
relations. To see this, we calculate the single-particle 
excitation spectrum, which is obtained as 

J 4(q,w) = J 4-(q,-w) + i4+(q,w) (4) 

with the photoemission (PES) spectrum 

A-(q,w) = £ K^- 1 !^!^}! 2 ^ - E^- 1 + <*)(5) 

va 

and inverse photoemission (IPES) spectrum 

A+(q,c) = J2 \{^ +1 \cU< h )\ 2 ^ - E^ +1 + <*)(6) 

where c qtT (c qCT ) is the Fourier transform of the creation 

(annihilation) operator c\ a (c, CT ). For the same reason as 
the calculation of the dynamical density-density correla- 
tion function, we employ the Lanczos method with the 
PBC where the operator is defined by 

Cqcr = — = e q z Ci a . (7) 
V i 

The calculation is carried out with a 18-site (L a = 6, 
L c = 3) cluster. 

4- Optical conductivity 

Finally, for more elaborate study on the behavior of the 
charge degrees of freedom and direct comparison with ex- 
periments, we calculate the optical conductivity defined 
by 

a ^ = ^Y / UJ \^ £} a \*l> }\ 2 6(u-E„ + E ), (8) 



with a component of the dipole operator 

D a = -e^Ti ■ v a (ni - -), (9) 

i 

where v a is a unit vector parallel to the a-direction. It 
is known that the OBC is more feasible for calculation 
of the optical conductivity spectra because finite-size ef- 
fects have much smaller influence on the results in com- 
parison with a case using the PBCM- Hence, we apply 
the OBC for both the a- and c-directions and study a 
finite-size cluster L a = 8, L c = 6 with the dynamical 
DMRG (DDMRG) method. The DDMRG method is an 
extension of the standard DMRG method for the calcu- 
lation of dynamical properties. In the DDMRG calcu- 
lation, a required CPU time increases rapidly with the 
number of the density-matrix eigenstates m, so that we 
try to keep it as few as possible. Because the DMRG 
approach is based on a variational principle, we have to 
prepare a 'good trial function' of the ground state with 
the density-matrix eigenstates as much as possible. We 
therefore keep m ~ 600 to obtain the true ground state 
in the first five DDMRG sweeps and keep m ~ 350 to 
calculate the spectrum of the system. As a result, the 
maximum truncation error, i.e., the discarded weight, is 
about 1 x 10 -3 , while the maximum error in the ground- 
state is about AE /N - 2 x 10" 2 . 



III. RESULTS 
A. Ground-state phase diagram 

So far, it has been suggested that the system defined 
by Eq.((T]) has three types of CO phases, depending on 
the values of the intcrsite Coulomb repulsions. The CO 
patterns for each phase are shown in the inset of Fig. [5] 
(a) vertical stripe for V c <C V p , (b) diagonal stripe for 
V c V p , and (c) three- fold alignment for V c ~ V p ^> 
t p . Other likely patterns are also shown: (d) horizontal 
and (e) randomly aligned stripes, which are not realized 
in the ground state. A random alignment CO pattern 
is formed with any combination of alternately-occupied 
chains parallel to the c-axis. Hence, the diagonal CO 
pattern is a special case of the random alignment ones. 

First of all, we note that the CO is observed as a state 
with a broken translational symmetry in our calculations: 
actually, there are two or more degenerate ground states 
and one of the degenerate states is picked out as the 
ground state by an initial condition of the DMRG calcu- 
lation as we here apply the OBC along the i-direction. 
Thus, a CO state includes more than two kinds of sites 
with different hole densities, i.e., hole-rich and hole-poor 
sites. In Fig. [21 we show DMRG results of the hole den- 
sity (rijgrich) and double occupancy (digrich) at the hole- 
rich sites i G rich for (a) V c = 1, (b) V p = 1, and (c) 
V p = V c - When both V p and V c are small, the analysis 
shows that the hole density is uniform over the system, 
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diagonal CO 



FIG. 2: (Color online) Hole density (m) and double occu- 
pancy (dj) at the hole-rich sites for (a) V p = 1, (b) V c = 1, 
and (c) Vj, = V" c . The inset shows the schematic represen- 
tation of (a) the diagonal stripe, (b) vertical stripe, and (c) 
three-fold CO patterns. The solid circle denotes a site where 
the hole predominantly resides. The other configurations are 
also shown: (d) horizontal and (e) an example of randomly 
aligned CO patterns. 
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FIG. 3: (Color online) Ground-state phase diagram of the ex- 
tended Hubbard model defined on the anisotropic triangular 
lattice. We assume the values t c = and U = 10. The result 
is obtained from the DMRG calculations for the charge distri- 
butions with a broken translational symmetry. The diameter 
of the circles is taken to be proportional to the hole density 
at the hole-rich sites. 



i.e., (rii) = 0.5 for all sites i, and the system is in the 
metallic state. As V p and/or V c increase, the charge fluc- 
tuations toward the CO states are enhanced and we find 
that (nigrich) starts to increase at the CO phase bound- 
ary. In the CO phase, the charge fluctuations are weak- 
ened again as the CO state is stabilized, where (riig r i C h) 
approaches a value in the atomic limit t p — ► 0. 

Let us now evaluate the phase boundaries between the 
metallic and CO phases. When we increase V p from zero 
with fixing V c = 1 [see Fig. Ufa)], (n ieric h) begins to in- 
crease around V p = 2.5; accordingly, a critical strength of 
the vertical CO state is obtained as (V p , V£) ~ (2.5, 1). 
At V p > V p c , (ni 6r i C h) increases rapidly to unity but the 
double occupancy remains almost zero. This means that 
(nearly) half-filled chains are alternated with (nearly) 
empty chains and the charge fluctuations are totally sup- 
pressed in the entire region of vertical CO phase. In 
analogy with the case of the vertical CO state, we can 
estimate a critical strength of the diagonal CO state as 
(Vp,V c c ) ~ (1,2.75) [see Fig. Hfb)]. We should how- 
ever note that (njgRich) seems to approach unity more 
slowly than that in the vertical CO phase. This implies 
that the charge fluctuations in the diagonal CO state are 
rather stronger than those in the vertical CO state. The 
reason of this strong charge fluctuation in the diagonal 
CO phase is the existence of nearly energetically degen- 
erate states, i.e., the horizontal and an infinite number 
of the randomly aligned CO patterns, which are shown 
in Fig. [2jd) and (e), respectively. This situation has 
also been confirmed by the mean-field approximation,? 3 
variational Monte Carlo (VMC) method^ and strong- 
coupling study of the spinlcss models The energies of 



those CO patterns including the diagonal one are exactly 
equal to each other in the atomic limit t p — > 0. 

As for the three-fold CO phase, the situation is some- 
what different from the other CO phases. If a relation 
V c = V p is kept, the three- fold CO instability begins 
to appear around a critical point {V p ,V^) ~ (3,3) [see 
Fig. I2I C )]- Above the critical point, we find the rapid 
increase of not only the hole density but also the dou- 
ble occupancy at the hole-rich sites. This implies that 
the competition between the effects of U and V's is es- 
sential for the presence of the three-fold CO phase. We 
also note that the three-fold CO state is metallic in the 
sense that the Drude weight is nonzero^ The presence 
of metallic CO phase has previously been suggested for 
the extended square-lattice Hubbard models; simply, the 
pocket-like small Fermi surface appears by doping the 
long-range CO phases with holes or electrons without 
destroying the long-range CO.— Note that the ground 
state has both the onsitc hole-pairing and CO, which may 
be essentially the same as the coexistence of the s-wave 
superconducting and CO states in the two-dimensional 
negative-?/ Hubbard models 

In Fig. [31 we determine the ground-state phase dia- 
gram from the results of the hole density, where the phase 
boundaries arc drawn with dashed curves. We can make 
it certain that the three substantial CO phases are distin- 
guished from a metallic phase with uniform charge distri- 
bution. This phase diagram is basically consistent with 
that from the VMC calculation.— The main difference is 
that a metallic region is sandwiched between the verti- 
cal and three-fold CO phases in our results. On the other 
hand, the diagonal and three-fold CO phases are contigu- 
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FIG. 4: (Color online) Transition between diagonal and three- 
fold CO states. The solid circle denotes a site where hole 
predominantly resides. The arrows denote the flow of charges 
toward the other CO state. 



ous, in agreement with the VMC result. The appearance 
of the metallic region in-between may be explained as 
follows. The diagonal and three-fold CO phases may be 
contiguous because the transition can be derived by the 
flow of charges only via the nearest-neighbor hopping in- 
tegrals t p (sec Fig. S]), whereas a drastic charge redistri- 
bution is required in the transition from the vertical CO 
pattern to the three-fold CO pattern and therefore the 
ground state needs to go through the metallic one. We 
may also suggest that this 'charge-flowing' transition be- 
tween the diagonal and three-fold CO phases occurs in 
the presence of the three-fold (diagonal) charge fluctua- 
tions in the diagonal (three-fold) CO phase as is evident 
in the calculated dynamical density-density correlation 
functions (see the following subsection). We also note 
that the three-fold CO phase is shrunk with increasing 
U. 



B. Dynamical density-density correlation 
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FIG. 5: (Color online) Charge excitation spectra calculated 
for the 18-site cluster with the Lanczos exact-diagonalization 
method. Only the results of (it, 0) are obtained for the 16-site 
cluster. Spectra at (V P ,V ) = (1,5), (5,1) [left panel], (0,0), 
(3,3), and (5,5) [right panel] are shown. 



The dynamical density-density correlation functions 
calculated by the Lanczos method for several sets of V p 
and V c values are shown in Fig. [SJ Each parameter set 
corresponds to a different phase as follows: (V P ,V C ) = 
(1,5) to the diagonal CO phase, (5,1) to the vertical 
CO phase, (5, 5) to the three-fold CO phase, and (0, 0) 
to the uniform metallic phase. In addition, we plot the 
results for (V P ,V C ) = (3,3), where the system is in the 
non-CO metallic phase but with strong three-fold CO 
fluctuations. When (V P ,V C ) = (0,0), the model flU is 
equivalent to the square-lattice Hubbard model at quar- 
ter filling and no particular modes of charge fluctuations 
are developed; we therefore find the broad spectral fea- 
tures for all the momenta, which come basically from the 
particle-hole transition in the corresponding noninteract- 
ing system. With increasing the intersite interactions, 
however, the broad features change to the low-energy 
sharp peaks, indicating that the excitations turn to be 
collective-mode like. 

Generally, the enhancement of low-lying peak at a par- 
ticular momentum must be expected for the CO insta- 
bility. When the system is unstable to the vertical [di- 



agonal] CO phase, the particular momentum is located 
at q = (tt, 7t) [q = (tt, 0)]. Looking at the spectra for 
iy P i V c ) = (5, 1), a sharp peak appears around u> = at 
q = (7T, 7r) and there are nearly no peaks at the other 
momenta. This means that the system is unstable to the 
vertical CO state with (almost) complete charge dispro- 
portionation. This result is consistent with the results of 
the hole density given in Sec. Ill A. On the other hand, 
in the spectra at (V P ,V C ) = (1,5), we see not only a 
low-energy peak at (tt, 0) but also relatively high-energy 
peaks at the other momenta; most of the high-energy 
peaks concentrate at momenta with q x q y > and few 
peaks appear at momenta with q x q y < 0. By the axes 
rotation, we find the charge fluctuations mostly occur 
along the c-axis, which are associated with the nearly de- 
generate states, i.e., horizontal and randomly aligned CO 
phases. It is moreover striking that the intensity around 
q = (2-7r/3, 2ir/3) is significantly large, which means that 
the three-fold charge fluctuations are rather strong al- 
though the system is in the diagonal CO phase. 

Let us now discuss the spectral features for the three- 
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fold CO metallic state, of which the particular momen- 
tum is located at q = (2n/3, 2tt/3). The charge excita- 
tion spectra in the CO metallic state are currently not 
well understood, so that a good opportunity for studying 
them may be offered here. We thus investigate a process 
of gradual change in the spectra with increasing the inter- 
site Coulomb interactions Vs. Since the three-fold CO 
phase lies around V p ~ V c > 3, we increase V p (V c ) from 
to 5 with keeping a condition V p = V c . The spectra for 
(^P)K) = (0,0) have the broad spectral features for all 
the momenta as mentioned above and these features ba- 
sically remain unchanged as far as V p and V c are less than 
~ 2. Only around (V^,, V c ) = (3, 3), we can recognize the 
changes clearly in the calculated spectra: the low-energy 
spectrum at q = (2tt/3,2tt/3) is particularly enhanced 
as anticipated, but those at q = (tt, 0) and q = (tt, tt) are 
also enlarged, while the spectral weights for the other 
momenta are strongly suppressed. We understand that 
these changes suggest that the three types of CO insta- 
bility are simultaneously developed and that, considering 
that the energies of the lowest-lying excitations for the 
three momenta arc nearly equal, the charge fluctuations 
for the three CO patterns are strongly competing. 

After increasing the intersite Coulomb interactions to 
(V P ,V C ) — (5,5), even greater changes are seen in the 
spectral features. At q = (2tt/3, 2tt/3), we find the en- 
hancement of a sharp peak around lu = 0, which indicates 
the strong charge fluctuations of the three-fold periodic- 
ity and is consistent with the results of the hole density. 
It is particularly worth noting that the low-energy spec- 
tral intensities at (it, 0) is continuously enhanced [and 
those at (tt, it) are suppressed]. 

Let us discuss some experimental situations here. We 
may point out first that the wave vectors qi and q2 
of the two diffuse peaks observed experimentall y 14 ^ 5 in 
the CO phase of 0-(ET)MM'(SCN) 4 [M =Cs, M =Co 
or Zn] should be equivalent to our momenta (tt, 0) and 
(27r/3, 27r/3), respectively. We may therefore suggest 
that the appearance of the two enhanced peaks in our cal- 
culations correspond to the coexistence of the two differ- 
ent charge fluctuations. Note that these enhanced peaks 
are already obtained for more realistic parameter values 
(V P ,V C ) = (3,3) in our calculations. It is also striking 
that the low-energy peaks appear in the entire Brillouin 
zone and the charge fluctuations are induced in all the 
momentum transfers. This situation is quite different 
from those for the diagonal and vertical CO stripes where 
the low-energy peak appears only at the particular mo- 
menta. 



C. Single- particle excitation spectrum 

In the diagonal and vertical CO phases, we just find 
the insulating band structures reflecting the CO states 
in our calculations; namely, the PES and IPES peaks 
are separated by the charge gap and the Fermi level lies 
between the two spectra. We thus focus on the spec- 
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FIG. 6: (Color online) Single-particle excitation spectra 
calculated for the 18-site cluster with the Lanczos exact- 
diagonalization method. The dashed lines denote the Fermi 
levels and the thin lines mark the dispersion of bands. 



tral features in the three-fold CO phases here. In Fig. O 
we show the single-particle excitation spectra calculated 
with the Lanczos method for (V p , V c ) = (0, 0), (3, 3), and 
(5,5). At (V p , V c ) = (0,0), we can clearly see the spec- 
tral weight separated into two bands due to the onsite 
Coulomb interaction, i.e., lower Hubbard band (LHB) 
and upper Hubbard band (UHB) with the distance ~ 10 
(= U). The LHB, which corresponds to the conduction 
band, can be approximately expressed by the dispersion 
of nonintcracting holes, 



e q = — 2t p cos q x — 2t p cos q y . 



(10) 



On the other hand, the UHB has relatively small spectral 
weight and narrow band width, which are general fea- 
tures in the strong-coupling regime. At (V p , V c ) = (3, 3), 
the shape of conduction band remains almost unchanged 
from that at (V p , V c ) = (0, 0). However, we can see that 
an overlap between the LHB and UHB occurs due to 
the competition between the onsite and nearest-neighbor 
Coulomb interactions. As a result, the band structure is 
seen to be reduced to a single dispersion line. The inter- 
action between holes is effectively weakened (or vanished) 



7 



by the frustration of the Coulomb interactions, so that 
the dispersion relation looks as if it were for the single- 
orbital model or the nonintcracting case. 

Surprisingly, the spectral features are drastically 
changed at (V P ,V C ) = (5,5). The spectral weight ap- 
pears to be separated to form several bands again. It is 
particularly worth noting that the width of conduction 
band is much narrower than that at (V P ,V C ) — (3,3). 
In the three-fold CO phase, the holes can conduct only 
through the hole-rich sites but the hole-rich sites are not 
connected by direct hopping integrals; in other words, the 
carriers form a small Fermi surfaced Thus, the conduc- 
tion band is centered at uj ~ 10(= U) and the band width 
should be very narrow. It may be related to the appear- 
ance of low-energy peaks in the entire Brillouin zone ob- 
served in the dynamical density-density correlation func- 
tion as discussed above. Let us then consider the peaks in 
the higher-energy range, uj 15. These peaks are essen- 
tially formed in connection with the empty (hole-poor) 
sites. An empty site is surrounded by three hole- rich sites 
and a hole-rich site is occupied by one or two holes. When 
a hole is added on an empty site, the excitation energy 
is uj ~ 2V P + V c = 15 if the surrounding three hole-rich 
sites are singly occupied and it is uj ~ AV P + 2V C = 30 
if the surrounding three hole-rich sites are doubly occu- 
pied. This is the reason why there are continuum weights 
from uj ~ 15 to 30 in the spectra for all the momenta. In 
addition, we can find the dispersion with bandwidth ~ 4 
around uj = 20. This dispersion is formed from a empty 
zigzag chain along the a-axis, which is nearly separated 
from each other by the hole-rich sites. 



D. Optical conductivity 

The optical conductivity spectra calculated by the 
DMRG method for several sets of V p and V c are shown 
in Fig. [71 For all cases, the spectral intensity for the in- 
direction seems to be a few times larger than that for 
the c-direction, which just reflects the difference between 
the interatomic distances d a : d c = V3 : 1 in the a and c 
directions. This leads to the intensity ratio I a : I c = 3 : 1. 

The most easily comprehensible result should be the 
ones in the vertical CO phase, i.e., {V p ,V c ) = (4,1) 
[see Fig. UJb)]. The shapes of spectra for the a- and 
c-dircctions are essentially the same, although the in- 
tensities are different as mentioned above. For both 
of the spectra, the lowest excitation appears around 
uj ~ 2.5, which corresponds to the optical (insulating) 
gap, and the other excitations lie in much higher fre- 
quencies, uj > 4. Also, the total weight of the spectrum 
is smaller than that in the other CO phases. These re- 
sults are consistent with the fact that the vertical CO 
state is rather solid and the charge fluctuations are quite 
small. The optical gap in the vertical CO phase may be 
estimated as ~ 3V P — 2V C — W with the bare bandwidth 
W = 8. 

Let us then turn to the diagonal CO phase, i.e., 




FIG. 7: (Color online) Optical conductivity spectra calculated 
with the DDMRG method. a(= a, c) is the direction of the 
electric field defined in Fig. [l] 



(V P ,V C ) = (1,4) [see Fig. Ufa)]. We find that the result 
for the a-direction is apparently different from that for 
the c-direction. The spectrum for the a-direction consists 
of the first-excitation peak at uj ~ 2.2 and high-energy 
broad features, which is similar to that in the vertical 
CO phase. This result reflects the solid CO and weak 
charge fluctuations along the a-direction, as shown in 
the dynamical density-density correlation functions. On 
the other hand, most of the spectral features for the c- 
direction appear in the lower-energy range. The optical 
gap is of the order of ~ 0.5, which is much smaller than 
that of the vertical CO phase. This is so because the 
diagonal CO state can be rearranged to the randomly 
aligned (or horizontal) CO state with small excitation 
energy. Thus, the charge fluctuations along the c-axis 
are much stronger than those along the a-axis. Since the 
energy difference between the diagonal CO phase and 
random alignment CO state is hardly changed even if V c 
increases, the low-energy structure is expected to be seen 
for any V c value. 

Lastly, we discuss the spectral aspects in the metallic 
regime around the three- fold CO instability [see Fig.[7Jc)- 
(e)]. With the OBC, the optical conductivity at uj = is 
ruled out and the lowest-energy scale is limited to an or- 
der of AW/L a L b ~ 0.67.-^ Therefore, the Drude spectral 
weight appears around the minimum excitation energy. 
At (V p , V c ) — (2, 2), the spectra for the a- and c-directions 
are quite similar. For each spectrum, small peak struc- 
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ture around lj ~ 0.7—0.8 should correspond to the Drude 
peaks. The spectral weight exists continuously from the 
'Drude frequency' to higher excitations and reaches the 
maximum around ui = 2, and very little spectral weight 
exists at higher-energy range, u > 3. At (V p , V c ) = (3, 3), 
we can find some tendencies to the three-fold CO state; 
in the spectrum along the a-axis the intensity around 
u> = 2 is about half of that at (V P ,V C ) — (2,2) and the 
high-energy broad features are relatively large. A large 
peak around ui = 2 at (V p , V c ) = (2, 2) splits into two 
peaks that indicate the appearance of individual excita- 
tions associated with the COs. The spectrum along the 
a-axis remains almost unchanged, except the Drude-like 
peak which becomes much smaller. At (V P ,V C ) = (4,4), 
we can clearly see the characteristic features of the three- 
fold CO state: in the spectrum along the a-axis the 
weight is separated into the low-energy Drude-like peak 
at w ~ 0.7—0.8 and high-energy broad features. We 
can understand that the former comes from the small 
Fermi surface in the three-fold CO phase and the latter 
comes from the excitations due to the three-fold charge 
fluctuations. On the other hand, the weight of the Drude- 
like peak is nearly zero in the spectrum along the c-axis. 
This means that the CO is more solid along the c-axis 
than along the a-axis. We can explain this as follows: 
3/2-filled chains are alternated with zigzag empty chains 
along the c-axis [see the inset of Fig. [5] (c)], whereas each 
chain along the c-axis is quarter filled, which means that 
the holes distribute uniformly along the a-axis. As men- 
tioned above, the optical conductivity along the a- (c-) 
axis is concentrated in the lower-energy range when the 
three-fold (diagonal) charge fluctuations are enhanced. 
We suggest that these features in the optical conductiv- 
ity spectra should give a good criterion to examine the 
CO patterns experimentally. 



IV. CONCLUSIONS 

The extended Hubbard model at quarter filling defined 
on the anisotropic triangular lattice have been studied 
by the Lanczos exact-diagonalization and DMRG meth- 
ods. We determine the ground-state phase diagram based 
on the results of the hole density and double occupancy. 
In the phase diagram, there exist three substantial CO 
phases (diagonal, vertical, and three-fold) and a metal- 
lic phase with uniform charge distribution. We find 
that the charge fluctuations (or instability to CO) with 
the three-fold periodicity appear for realistic parameters 
V p ~ V c ~ 0.3Z7. 

We suggest that the transition between the diagonal 
and three-fold CO phases is derived by the flow of charges 
only via the nearest-neighbor hopping integrals, which 
may lead to the coexistence of the diagonal and three-fold 
charge fluctuations near the phase boundary. The coex- 
istence features arc also found in the density-density cor- 
relation functions; in the diagonal (three-fold) phase, not 
only the particular peak but also the low-energy peaks 



corresponding to the three-fold (diagonal) charge fluc- 
tuations are enhanced. In the vertical CO phase, only 
the low-energy peak at the particular momentum is vis- 
ible and there is no charge fluctuations associated with 
the other CO states. Moreover, we find that there ap- 
pear the low-energy collective-modes-like excitations in 
the entire Brillouin zone when the three-fold charge fluc- 
tuations are very strong. The wave numbers of the diag- 
onal and three-fold charge fluctuations arc equivalent to 
those of the two X-ray diffuse peaks qi = (0, k, 1/2) and 
q 2 = (2/3, k, 1/3) measured in 6»-(ET)CsM'(SCN) 4 with 
M =Co or Zn. If we could make the dihedral angle larger, 
the value of V p /V c decreases and only the diagonal CO 
fluctuation remains. This is consistent with the fact that 
only a diffuse peak with wave number qi is observed in 
the compounds with M =T1 and Rb which have a larger 
dihedral angle than those with M =Co or Zn do. 

We also study the single-particle excitation spectrum 
to see the evolution of the fundamental low-lying excita- 
tions in the presence of strongly frustrated correlations. 
When the onsite and intersite Coulomb interactions com- 
pete (V p ~ V c ~ U/3), the interaction between holes is 
effectively diminished and the dispersion relation looks 
as if it is of the noninteracting case. When the two inter- 
site Coulomb interactions compete (V p ~ V c > U/3), the 
carriers form very narrow conduction band with a small 
Fermi surface. 

Furthermore, we examine the optical conductivity 
which reflects the characteristic features for each CO 
phase. In the vertical CO phase, the spectra for both the 
directions explicitly represent the insulating features and 
the optical gaps are large, estimated as ~ 3V p — 2V C — 8. In 
the diagonal CO phase, the spectrum for the a-direction 
is essentially the same as that in the vertical CO phase 
but that for the c-direction is located in the lower-energy 
range. This is so because the diagonal CO state can 
be rearranged to the randomly aligned (or horizontal) 
CO state with small excitation energy. In the three- 
fold CO phase, the spectra indicate the presence of the 
separated low-energy Drude-like peak and high-energy 
broad features, leading to the anomalous metallic states 
in the system. The optical conductivity along the a- (c-) 
axis is concentrated in the lower-energy range when the 
three-fold (diagonal) charge fluctuations. These spectral 
features may give a good criterion to examine the CO 
patterns experimentally. Finally, we make a comment 
on the difference in the experimental spectra between 
MM' =CsZn and MM' =RbZn. Because we find that 
an increase in V p /V c makes the diagonal CO state unsta- 
ble, it is the smallness of t c that is essential for the strong 
insulating features in the spectra of MM' =CsZn. 

We hope that our studies will help us understand the 
mechanism of charge fluctuations in the extended Hub- 
bard model defined on the anisotropic triangular lattices 
and hence will offer useful suggestions to some aspects of 
the three-fold charge fluctuations observed in 0-(ET)2X. 
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